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ABSTRACT 

The feedback between ice-snow albedo and temperature is 
included in a one-dimensional radiative-convective climate 
model. The effect of this feedback on sensitivity to changes 
in solar constant is studied for the current values of the solar 
constant and cloud characteristics. The ice-snow albedo feed- 
back amplifies global climate sensitivity by 33% and 50%, respec- 
tively, for assumptions of constant cloud altitude and constant 
cloud temperature. 

Introduction. One-dimensional radiative-convective models 
which determine the thermal structure of the atmosphere as a result of 
the balatioe between the radiative flux and a parameterized convective 
flux are useful tools for climate studies. Such models, which can include 
realistic vertical distributions of radiatively-important atmospheric con- 
stituents, can be used to examine the role that these constituents play in 
determining the global mean temperature structure. One major draw- 
back of these models is their lack of the ice-snow albedo-temperature 
feedback which is of great importance because of the large difference in 
albedo between ice and ice-free areas. The present work describes a 
simple method of incorporating the albedo-temperature feedback in a 
one-dimensional radiative-convective climate model. In addition the ef- 
fect of this feedback on the sensitivity of the climate to changes in solar 
constant is studied and compared with the sensitivity found in other mod- 

Analysis. The one-dimensional radiative-convective model des- 
cribed in Wang and Lacis (1979) is used in the present study. Jn order to 
incorporate the ice-snow albedo-temperature feedback in a global model, 
it is necessary to know the response of the global mean surface albedo to 
changes in the ice-snow line as a result of changes in global mean sur- 
face temperature. 

We assume first the zonal mean surface albedo as l^^s a simple 

a 0 < X < Xg 
b x<x<l 

fcj * . -- — 

where a and b are constants characterizing the ice-free and ice-covered 
surface albedos, respectively, and xs is the sine of the latitude of the 
edge of the ice sheet. For the ice-sheet edge, the prescription by Budyko 
(1969) is adopted, i. e. , if the zonal mean surface temperature Ts is less 
than -10° C, ice will be present. We will further use the simple formu- 
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latlon of the latitudinal temperature structure used by North (1975). In 
this formulation the temperature structure is approximated by the first 
two Legendre pol 3 momials in a series expansion, i.e. , 




( 2 ) 


2 

where P 2 (^) = (3x -l)/2. The first term T is the global mean surface 
temperamre since the second term vanish^ upon integration from 0 to 1. 
This two-mode approximation with T|2«-28.40 C yields a reasonably good 
temperature distribution compared with observations. 

In the atmosphere To varies with the tropospheric lapse rate be- 
cause of the effect of large scale eddies (Stone, 1978). Since the tropo- 
spheric lapse rate is essentially constant in the radiative-convective mo- 
del, we can adopt the simple parameterization suggested by Stone (1978) 
for taking into account the effect of the eddies on temperature structure 
in climate calculations, i.e. , To = constant = -28.4 Q Then setting 
x = Xg, T = Tg = -10° C In Equation (2) we obtain 

Xg = (0.5681 + 0.02347 T . (3) 


This equation describes how the extent of the polar ice cap changes as 
the global mean surface temperature changes. 

Next, we wish to determine the relationship between the global 
mean surface albedo and Xg. Employing Equation (1), we obtain 


•'n 


Sg(X, 


X ) 

G' 


Og(X,Xg)dX 




Sg(x,Xs)dx 


( 4 ) 


where Sg is the zonal mean solar radiation received at the surface. Ss 
is a complicated function, which depends on solar zenith angle and at- 
mospheric transmissivity. As an initial calculation, to test the method, 
we assume that Sg has the same latitudinal distribution as the solar inso- 
lation at the top of the atmosphere, Sq. This approximation will be relaxed 
ultimately. Also we will use North ^s (1975) two mode approximation for 
Sq, i.e., we choose 

Sg(x,Xg) = [1+S2P2 (x)]S (5) 


where S is the mean value of Sg and Sg = -0.482, Substituting Equations 
(1) and (5) into (4), we obtain the following expression 

«g = [b(l-Xg) + axg ] + S2(b-a)(Xg-Xg )/2 (6) 


where the first term is the geometrical albedo. Equations (3) and (6) 
when coupled give a highly non-linear relation describing the feedback be- 
tween ice-snow albedo and temperature. 

Results and discussion. First we examine the sensitivity of the 
present radiative-convective model with and without the surface albedo- 
temperature feedback. A convenient measure of the sensitivity of the 
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Table !• Comparison of the model sensitivity parameter g and the ice- 
snow albedo amplification parameter y for several climate %iodels. 
is the dependence of zonal surface-atmosphere albedo upon zonal mean^ 
surface temperature (cf. , Sellers, 1969). 


MODEL 


y 

A. Energy balance models 

Budyko (1967) 

155 

1.58 

Lianand Cess (1977) 

147 


C^(x) 


0. 25 

- 0.009 


1.68 

= 0.004 


0. 27 

Sellers (1969) 

150 


= 0.009 


1,17 

Gal-chen and Schneider (1976) 

134 


= 0.009 


0.98 

= 0.004 


0.22 

B. Radiative-convective models 

Wetherald and Manabe (1975) 

I 

128 


II* 

114 


Present model 

FCAt 

112 

0.33 

FCT + 

140 

0.49 

♦Uses the radiation scheme of Rodgers and Walshaw (1966). 


+ FCA denotes fixed cloud altitude. FCT denotes fixed cloud temperature. 



Solar Constont, 5/5. 

Fig. 1. Position of the edge of the polar ice cap as a function of the so- 
lar constant, as determined by the present model (solid curve) and by 
North (dashed curve). 
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( 7 ) 


global climate model is the parameter defined as 

6 = 




dS 


The enhancement in global sensitivity due to surface albedo-temperature 
feedback can be described by the amplification parameter, 




( 8 ) 


where B is the model sensitivity in the absence of the surface albedo- 
temperaftire feedback. Values of Bq and Y for the present model and 
several other climate models are given in Table 1. In the calculations 
the present value of the solar constant S was chosen to be 1365 W m"^ 
and values of a and b were chosen to be 0. 09 and 0. 6, respectively. The 
characteristics of the model corresponding to the present climate are: 
a g = 0.1, Xg = 0.964, Tg = 15.37° C and global mean planetary albedo 
0.306. 

In general, radiative-convective models yield smaller values of 
B than energy balance models, since the former do not include the 
ot - T positive feedback. The most accurate values of B Y 
tfioseigiven by Lian and Cess (1977) for whose calculations 

Included the most realistic parameterizatlons ofxhls feedback. The re- 
sults of our model with this feedback included are also listed in Table 1. 
The enhancement of the model’s sensitivity is in good agreement with 
Lian and Cess’ results, although the magnitude of the enhancement is 
subject to modification when a more accurate representation of Sgls used. 

We also calculated the variation of the ice line with solar constant 
in the radiative-convective model. The solution is Illustrated in Fig. 1. 
The qualitative behavior Is similar to that in the 1-D heat balance models. 
As the solar constant is decreased, a point is reached at which the ice 
cap is unstable and will grow spontaneously until it covers the whole 
earth. For the approximation to Sg we have used, this occurs when the 
solar constant has decreased by 12%. 

We conclude that the surfcice albedo-temperature feedback can be 
included in a one-dimensional radiative -convective climate model in a 
qualitatively realistic way. This addition represents a significant im- 
provement in the ability of such models to simulate the effects of radia- 
tive perturbations on climate in a realistic way. 
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